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Q Abstract — In wireless ad hoc networks, distributed nodes 
can collaboratively form an antenna array for long-distance 
communications to achieve high energy efficiency. In recent work, 
Ochiai, et al., have shown that such collaborative beamforming 
can achieve a statistically nice beampattern with a narrow main 
lobe and low side lobes. However, the process of collaboration 
introduces significant delay, since all collaborating nodes need 
access to the same information. In this paper, a technique that 
significantly reduces the collaboration overhead is proposed. It 
consists of two phases. In the first phase, nodes transmit locally 
in a random access fashion. Collisions, when they occur, are 
viewed as linear mixtures of the collided packets. In the second 
phase, a set of cooperating nodes acts as a distributed antenna 
system and beamform the received analog waveform to one or 
more faraway destinations. This step requires multiplication of 
the received analog waveform by a complex number, which is 
independently computed by each cooperating node, and which 
enables separation of the collided packets based on their final 
destination. The scheme requires that each node has global 
knowledge of the network coordinates. The proposed scheme can 
achieve high throughput, which in certain cases exceeds one. 

Keywords: distributed wireless systems, cooperation, 
beamforming 



I. Introduction 

Energy is often a scarce commodity in wireless sensor 
networks, as wireless sensors typically operate on batteries, 
which in many cases are hard to replace. Similarly, due to 
cost consideration, nodes in wireless ad hoc sensor networks 
are commonly equipped with only a single omnidirectional 
antenna. Thus, in order to transmit information over long 
distances while conserving energy and maintaining a certain 
transmission power threshold, multihop networks have been 
the preferred solution. However, there are several challenges in 
transmitting real-time services over multiple hops. For exam- 
ple, the traditional CSMA/CA based media access control for 
avoiding collisions does not work well in a multihop scenario 
because transmitters are often out of reach of other users' 
sensing range. Thus, as packets travel across the network, 
they experience interference and a large number of collisions, 
which introduces delays. Also, multihop networks require a 
high node density which makes routing more difficult and 
affects the reliability of links [1]. 

1 This work was supported by the National Science Foundation under Grants 
ANI-03-38807, CNS-06-25637 and CNS-04-35052, and by the Office of Naval 
Research under Grant ONR-N-00014-07-1-0500 



Recently, a collaborative beamforming technique was pro- 
posed in [3], in which randomly distributed nodes in a network 
cluster form an antenna array and beamform data to a faraway 
destination without each node exceeding its power constraint. 
The destination receives data with high signal power. Beam- 
forming with antenna arrays is a well studied technology; 
it provides space-division multiple access (SDMA) which 
enables significant communication rate increase. The main 
challenge with implementing beamforming based on randomly 
distributed nodes in the network is that the geometry of 
the network changes dynamically. It was shown in [3] that 
such a distributed antenna array can achieve a nice average 
beampattern with a narrow main lobe and low side lobes. 
The directivity of the pattern increases as the number of 
collaborating nodes increases. Such an approach, when applied 
in the context of a multihop network reduces the number of 
hops needed, thereby reducing packet delays and improving 
throughput. However, one must take into account the overhead 
required for node collaboration, i.e., sharing of the information 
to be transmitted jointly by the nodes. If a time-division 
multiple-access (TDMA) scheme were to be employed, the 
information-sharing time would increase proportionally to the 
number of nodes involved in the collaboration. 

In this paper we present a technique based on the idea 
of collaborative beamforming, that reduces the time required 
for information sharing. The technique also allows different 
nodes in the network to transmit simultaneously. Collaborating 
nodes receive linear mixtures of the transmitted packets. 
Subsequently, each collaborating node transmits a weighted 
version of its received signal. The weights are such that one or 
multiple beams are formed, each focusing on one destination 
node, and reinforcing the signal intended for a particular des- 
tination as compared to the other signals. Each collaborating 
node computes its weight based on the channel coefficients 
between sources and itself, estimated by orthogonal node 
IDs embedded in the packets. The proposed scheme achieves 
higher throughput and lower delay with the cost of lower SINR 
as compared to [3]. 

II. Background on collaborative beamforming 

For simplicity, let us assume that the sources and destina- 
tions are coplanar. We index source nodes using a subscript i, 
with ti denoting the z-th node. The locations of these nodes 
follow a uniform distribution over a disk of radius R, and it is 



assumed that each node knows its own location. We denote the 
location of t. L in polar coordinates with respect to the center 
of the disk by (ri,ipi). 

Suppose that a set of N nodes designated as collaborating 
nodes have access to the same signal, i.e., s(n), whose desti- 
nation is at azimuthal angle 4>q. Let d,(0o) denote the distance 
between the z-th collaborating node and the destination. The 
initial phases at the collaborating nodes are set to 

*i(0o) = -ydiOfo), i=l,...,N (1) 

which requires knowledge of distances relative to wavelength 
between nodes and destination, and applies to the closed-loop 
case [3]. Alternatively, the initial phase of node i could be: 

*i(<M = -r-n cos(>o - ipi) (2) 

which requires knowledge of the node's position relative to 
some common reference point, and corresponds to the open- 
loop case [3]. In both cases synchronization is needed, which 
can be achieved via the use of the Global Positioning System 
(GPS). 

The channels between collaborating nodes and destination 
are assumed to be idential for all nodes. The corresponding 
array factor given the collaborating nodes at radial coordinates 
r = [ri,...,nv] and azimuthal coordinates ip = V'Af] 
is 
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Under far-field assumptions, the array factor becomes [3]: 
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{ri/R) sin(V>i - 
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\{4>q + 0)) with the following pdf: 



-1< z < 1 



Finally, the average beampattern can be expressed as [3] 

PavWO = E Z {\F(^)\ 2 } 
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where Ji(.) is the first-order Bessel function of the first kind. 
When plotted as a function of <f>, P w {4>) exhibits a main lobe 
around </>o> and side lobes away from 4>o- It equals one in the 
target direction, and the sidelobe level approaches 1/N as the 
angle moves away from the target direction. The statistical 
properties of the beampattern were analyzed in [3], where it 
was shown that under ideal channel and system assumptions, 
directivity of order N can be achieved asymptotically with N 
sparsely distributed sensors. 

As we have noted, all of the collaborating nodes must have 
the same information to implement beamforming. Thus, the 
active nodes need to share their information symbols with all 



collaborating nodes in advance. If a time-division multiple- 
access (TDMA) scheme were to be employed, the information- 
sharing time would increase proportionally to the number of 
active nodes. In the following, we propose a novel scheme to 
reduce the information-sharing time and also allow nodes in 
the network to transmit simultaneously. 

III. The proposed scheme 

Here we refine the model of [3], focusing more directly on 
the physical models for the signal, fading channel and noise. 
Besides the assumptions in Section [TT] we will further assume 
the followings: 

(1) The network is divided into clusters, so that nodes in 
a cluster can hear each other. In each cluster there is a node 
designated as the cluster-head (CH). Nodes in a cluster do not 
need to transmit their packets through the CH. 

(2) A slotted packet system is considered, in which each 
packet requires one slot for its transmission. Perfect synchro- 
nization is assumed between nodes in the same cluster. 

(3) Nodes operate under half-duplex mode, i.e., they cannot 
receive while they are transmitting. 

(4) Nodes transmit packets consisting of PSK symbols hav- 
ing the same variance a 2 . Each transmitted packet contains (in 
fixed locations) a set of pilots comprising the user ID, followed 
by a set of pilots comprising the destination information. 

(5) Communication takes place over flat fading channels. 
The channel gain during slot n between source and collab- 
orating node Cj is denoted by ay (n). It does not change within 
one slot, but can change between slots. The gains a,j follows a 
Rayleigh fading model, being i.i.d. complex Gaussian random 
variables with zero means and variances a 2 . 

(4) ^ ^ ne com Pl ex baseband-equivalent channel gain between 
nodes Ci,dj is [4]: = b ■ ei~^~ where is the distance 
between nodes Ci,dj, and where b is the path loss. The 
distances between collaborating nodes and destinations are 
much greater than distances between source and collaborating 
nodes. Thus, b is assumed to be identical for all collaborating 
nodes and equals the path loss between the origin of the disk 
over which the nodes are distributed to the destination. 

(7) the noise vectors are uncorrected, complex, zero-mean 
white Gaussian vectors. 

Suppose that cluster C contain J nodes. At slot n, nodes 
ti,. .. ,tn need to communicate with node d\ , . . . , dn that be- 
long to cluster C\, . . . , Ck, respectively. The azimuthal angle 
of destination dj is denoted by <pi . The packet transmitted by tj 
consists of M symbols &j(n) = [sj(n; 0), . . . , Sj(n; M — 1)]. 
Due to the broadcast nature of the wireless channel, non-active 
nodes in cluster C hear a collision, i.e., a linear combination 
of the transmitted symbols. More specifically, node c, hears 
the signal 



(3) 
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x i( n ) = y^aji(n)sj(n) +Wj(n) 



(7) 



where w,(u) = [wi(n; 0), . . . , Wi(n; M — 1)] is noise vector 
with variance crj at the receiving node a. 



Once the CH establishes that there has been a transmission, 
it initiates a collaborative transmission period (CTP), by send- 
ing to all nodes a control bit, e.g., 1, via an error-free control 
channel. The CH will continue sending a 1 in the beginning 
of each subsequent slot until the CTP has been completed. 
The cluster nodes cannot transmit new packets until the CTP 
is over. 

Suppose that each transmitted packet includes an ID se- 
quence, so that IDs are orthogonal between different users. 
The channel coefficients can be estimated by cross correlating 
Xj(n) with known user IDs as in [2]. If the magnitude of 
the cross-correlation in greater than some threshold, then the 
corresponding user is in the mixture, and the value of the cross- 
correlation provides the corresponding channel coefficient. The 
information of destination nodes could be obtained in a similar 
way. 

Each node q uses cross-correlation operations with known 
orthogonal user IDs to determine which users are in the mix- 
ture, the corresponding destinations, and also the coefficients 

a,ij. 

Let d m denote the destination of s m (n). In slot n + m, m = 
1, . . . , K, each collaborating node c, transmits the signal: 

Xi(n + m) = x i {n)fj, m a mi (n)e~ : > 2 z~ di( -' l ' m) (8) 

where di(<p m ) denotes the distance between user i and des- 
tination node d m with azimuth <fi m , and /i m is a scalar to 
adjust transmit power and same for all collaborating nodes. In 
addition, ii m is of the order of 1/N. 

Given the collaborating nodes at radial coordinates r = 
[ri, rjv], azimuthal coordinates i/> = [?/>i, i/V] and the 
path loss b, the received signal at direction </> is: 

N 

y(<j>;m\r,ip) = Y^&x^n + m)e i ^ di ^ + v(n + m) (9) 

i=l 

where v(n + m) is the noise vector with variance a 2 at the 
receiver during slot n + m. 

Let us consider the received signal at the destination d m at 
angle <p m during slot n + m: 

N 

y(</> m ; m) = ^ 6x;(n + m) + v ro (n + m) 

N 

= b^2[\a mi (n)\ 2 s m (n) + a* nl (n)wi(n) 

K 

+ a *ni( n ) a 3i( n ) a j( n )] +v m (n + mXlO) 

3=1 

As N -> oo, jjEti\a ml (n)\ 2 -> E{\a ml \ 2 } = a 2 a . 
Also > ^Eti a ml (n)a 3l (n) -» E{a ml a 0l } = 0, due to 
the fact that for j ^ rn, the channel coefficients a m i, aji are 
uncorrelated and have zero mean. Finally, as N — > oo and 
omitting the noise, y(4> m ; m|r, tp) — > N ^ m b m <jls m (n). Thus, 
the destination node d m receives a scaled version of s m (n). 



The beamforming step is completed in JsT slots, reinforcing 
one source signal at a time. 

Assuming that all of the K source packets have distinct 
destinations at different resolvable directions, multiple beams 
can be formed in one slot, each beam focusing on one direction 
and reinforcing one source signal. The transmitted signal from 
each collaborating node would be: 

K 

x i (n + l)=Xi(n)^/ima^(n)e-^ d *^ (11) 

m=l 

The received signal at destinations would then be: 

yO|r, -0) = bSii(n + l)e j ^ d ' w + v(n + 1) (12) 

It can be shown that as N — > oo and omitting noise, 
y((f) m \r, if>) — > iV/^ m 6cr^s m (n), for m = 1, K. Thus, each 
of the K beams transmits a scaled version of a source signal 
to its destination. Based on (O, ( TTZb . we can easily extend the 
mathematical formulation to the scenario that not all of the K 
packets have distinct destinations. 

In the rest of the paper, for simplicity we will consider 
only the case in which a single beam is formed during slot 
n + rn, focusing on destination d m . The results obtained 
under this assumption can be readily extended to multiple 
simultaneous beams. The time index can be omitted without 
causing confusion. Also, we are particularly interested in the 
average properties over a statistical ensemble, so the analysis 
can be based on one sample of a packet assuming samples are 
independent. Substituting y, x, x, w, v in ©-(O by y, x, x, 
w, v (i.e., with one of their samples) respectively, we get: 

JV 

y(4>;m\r,il>) = ^baae? 3 ^^ +v (13) 
i=i 

IV. Performance of the average beampattern 

In this section we analyze the average beampattern. Under 
the far-field assumption, and following the steps in [3], ( fT3b 
can be expressed as: 

N K 

y{<j>; m\v, 0) = £ u m ba* ml e~^^ (]T a jiSj + w t ) + v 

i=l j=l 

(14) 

where a(<j>; </> m ) and Z{ are the same as those in (0]). 

The far-field beampattern or the received power is defined 

as: 

P(<j>) = \y(<j>;m\r 1 xP)\ 2 (15) 

Taking into account the assumptions on the channel coeffi- 
cients, it can be readily shown that the average beampattern 
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B. Transmit Power and Average SINR 

In this section, we analyze the SINR under the same 
transmit power for the scheme of [3] and the proposed scheme. 

Although the focus in [3] was on the statistical properties 
of the beampattern, we can extend ©-© to a physical model 
including signal, path loss and noise. The received signal of 
the destination d m is given by: 

N 

_ y(^;m\r,i,) = Y,bll m s m e-^^ +v (19) 

+ J2 i^i 2 ^i a ™i 2 i a ^i 2+ Zli a ™i 2 i Wi i 2 + h 2 } 1=1 

i =1 i— 1 i— 1 - 2tv j ii \ 

j* m ci simply transmits [i m s m e 3 >• a,Wm ', and the transmit 

Nn 2 m b 2 [2a 2 a* + (N - \)a 2 s a A a E{e- 3a{ ' i ''' (t ' m){z ^ z '- ) } power is thus p^a 2 . The average beampattern at the target 
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Then, 
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The term N n^b 2 a 2 a^(3 represents the average 
the sidelobes that is independent of the angle (j>. 
(TTTI i is of the similar form with (O. Thus, other 
of the average beampattern like peak/zero positions 
bandwidth/sidelobe region can be easily obtained 
corresponding results of [3], and are omitted here. 

V. Performance of the network 
A. Throughput 

Suppose that K packets need to be transmitted. For the 
collaborative beamforming scheme of [3], each packet must be 
shared by the N beamforming nodes. If one wished to avoid 
collisions, TDMA would be a natural way to implement infor- 
mation sharing. With the use of TDMA, in each slot one active 
node is scheduled to broadcast its packet to other nodes within 
the same cluster. The sharing of information would require K 
slots. Via the use of multiple beams, the beamforming to the 
destination would require 1 slot if destinations are distinct, or 
up to K slots in the worst case where all nodes have the same 
destination. Thus, for the scheme of [3], the throughput, T 
satisfies: 1/2 < T < K/(K + 1). 

In the proposed scheme, combinations of K packets enter 
the system and reach N collaborating nodes in 1 slot. Via 
multiple beams, the packets can be delivered to their destina- 
tions in one additional slot, if the destinations are distinct. The 
throughput is then K/(l + 1) = K/2. If two or more packets 
have a common destination, the beamforming will need to 
take more than one slot. In the worst case where all packets 
have the same destination, K slots will be needed, resulting 
in a total throughput of K/(l + K). Thus, K/(l + K) < 
^proposed < K/2. Note that the throughput of the proposed 
scheme could be greater than 1. 



= N 2 f m b 2 a 2 



2 \ and the SINR is 



SINR: 



(20) 



(21) 



For the proposed scheme, the collaborating node c s ; trans- 
mits Xi given by ( fTTT i. It can be shown the average power of 
Xi equals 



E{xiX*} = fJ, m <7 s <r a {3' 



(22) 



where /?' = K + 1 
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To keep the same average transmit power as ^„cr 

-2 



(23) 



Under this value of pb m and based on ( fTTI i. the average 
beampattern at <f> m is given by 



N 2 fl 2 m b 2 a 2 
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(24) 



Thus, the average received power of the proposed scheme 
(without the noise term a 2 ) is f3' times less than that of 
d20l >. In other words, each collaborating node in the proposed 
scheme needs to use j3' times more transmit power, in order 
to (asymptotically) achieve the same average received power 
as ©. 

From 10), (ED, (US, (ED, the SINR at the destination d m 
for the proposed scheme is 



SINR 
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(25) 



Thus, compared with (ED, the average SINR of the pro- 
posed scheme at the destination is still asymptotically scaled 
down by /3'. 



VI. Simulations 

A. simulation setup 

We investigate the performance of the above scheme for 
different numbers of collaborating nodes, i.e. N. For conve- 
nience, and without loss of generality, we divide the nodes in a 
cluster into two pools: one contains all potential active nodes, 
the total number of which is fixed J = 32 and another contains 
N collaborating nodes. The directions of destinations are 
uniformly distributed in [0, 2tt). The locations of collaborating 
nodes are uniformly distributed within a disk with R = 10 
(i.e., the radius normalized by the wavelength). The channels 
among nodes in a cluster are selected from zero-mean complex 
Gaussian processes, which are constant within one slot, but 
vary between slots. 

B. BER performance 

In the proposed scheme, noise enters the collaborating 
nodes, with variance cr^, and the destination, with variance 
<r^. Let us define 71 = CsCq/cJ, which represents the 
average SNR in the process of information sharing, and define 
72 = N 2 y?^ 2 a 2 s o\l 2 . Note that 72 is also independent of 
N since fi m is of the order of 1 /N. The overall SINR in 
can be rewritten by 



SINR 
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which is determined by 71, 72, K and TV. 

We fix 71 to 20 dB to investigate the performance of 
bit error rate (BER) for different values of 72 and N. We 
perform a Monte-Carlo experiment consisting of 10 6 repeated 
independent trials. Each packet contains BPSK symbols and 
K = 4 nodes are transmitting all the time. Fig. \V[a) shows 
the BER vs. 72 for the case in which only one beampattern is 
formed in each slot, for different values of N. Solid lines 
correspond to estimated channels, number of active nodes 
and destination information. One can see that BER decreases 
as 72 and N increase. Dashed lines correspond to perfect 
knowledge of channels, number of active nodes and destination 
information, and can be considered as lower bounds of BER 
performance. Fig. Ulb) shows the BER performance in which 
all of four simultaneous beams are always formed in one slot. 
Compared with Fig.QJa), more collaborating nodes are needed 
to achieve the same BER level. 

C. Throughput Performance 

To investigate performance under certain traffic load p, we 
perform a random experiment consisting of 1,000 repeated 
independent trials. In each trial, all users are statistically 
the same, and each user sends out packets with probability 
p/J. The throughput is defined as the average number of 
packets that were successfully transmitted in a time slot. Each 
packet contains 424 bits with QPSK symbols. The nodes' 
ID sequences are selected based on a Jth order Hadamard 
matrix and the IDs are used to determine the active nodes 
involved in collisions as well as to estimate the channels. A 



maximum likelihood decoder is used at the destination nodes 
to recover the symbols. Packets received at the destinations 
with BER higher than 0.02 are considered to be lost or 
corrupted. 71 = 72 — 20 dB. 

In Figs. |2] (a)-(c), we show the throughput performance 
for different values of p and N, allowing up to 1, 2, 3 
simultaneous beampatterns per slot, respectively. The curves 
with legend "ideal" correspond to the case where all the 
transmitted packets are successfully received, which can be 
considered to be an upper bound on the proposed scheme. 
One can see that the increase of N can result in throughput 
improvement. Furthermore, for small N, one should choose 
a small number of simultaneous beampatterns to improve 
throughput. Fig. |2|d) shows the throughput in which all the 
K beampatterns are allowed to be formed in one slot. Note 
that N = 1, 024 enables throughput of almost K/2. 

VII. Conclusions and Future work 

In this paper we have proposed a technique for reducing 
the time needed for information sharing during collaborative 
beamforming, and for allowing simultaneous transmissions. 
The proposed scheme can achieve high throughput at the cost 
of reduced SINR. An analysis for the average beampattern and 
network performance has also been provided. Our analysis is 
based on a number of ideal assumptions on the system. In 
future work, we plan to investigate the effects of imperfect 
channel/phase and non-identical path loss, and also to seek 
closed-form BER expressions. 
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Fig. 1. BER performance under different 72 and N (K = 4) 




(b) up to two beams per slot 
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Fig. 2. Throughput performance under different traffic load and N 



